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We present detailed derivations, various improvements and application to concrete experimental 
data of spin squeezing inequalities formulated recently by some of us [Phys. Rev. Lett. 95, 
120502 (2005)]. These inequalities generalize the concept of the spin squeezing parameter, and 
provide necessary and sufficient conditions for genuine 2-, or 3- qubit entanglement for symmetric 
states, and sufficient condition for general JV-qubit states. We apply our method to theoretical 
study of Dicke states, and, in particular, to VF-states of N qubits. Then, we analyze the recently 
experimentally generated 7- and 8- ion W-states [Nature 438, 643 (2005)]. We also present some 
novel details concerning this experiment. Finally, we improve criteria for detection of genuine 
tripartite entanglement based on entanglement witnesses. 



I. INTRODUCTION 

Experimental generation and characterization of en- 
tanglement on a macroscopic, or mesoscopic scales seem 
to be one of the necessary prerequisites of scalable quan- 
tum information processing. A spectacular progress has 
been achieved recently in the area of quantum correlated 
systems of atoms, and in particular macroscopic atomic 
ensembles [lj . The main goal of these studies is to achieve 
an efficient quantum interface between light and atoms 
with spin, or pseudo-spin internal states, using the gener- 
alized quantum Faraday effect. Such settings already al- 
lowed one to demonstrate entanglement of distant atomic 
objects 0, or deterministic memory for light Q that can 
be retrieved using quantum teleportation (4j. Entangle- 
ment between light and atoms, and between atoms them- 
selves plays, of course, essential role in these experiments. 

It worth stressing that the light-atoms interface based 
on using the quantum Faraday effect does not only allow 
one to measure and detect atomic states. It does also pro- 
vide a tool for manipulations and engineering of quantum 
fluctuations of atomic spins. The latter possibility might 
be of fundamental importance for the future implementa- 
tions of distributed quantum information processing. In 
particular, the methods of atomic ensembles can be car- 
ried over to another rapidly developing area of ultracold 
atomic gases. Here, the interest would be to measure 
characterize, and finally engineer quantum fluctuations 
of the total atomic spin in spinor ultracold gases (for a 
review see [f| ) that has been intensively studies since the 
seminal theory papers of T.-L. Ho @ and T. Ohmi and 
K. Machida [3], as well as the experiments performed 
by the MIT group on optically trapped sodium Bosc- 
Einstein condensates (BEC) Q. Particularly interesting 
are prospects of applications of these methods to strongly 
correlated states of spin ultracold gases in optical lattices 
1- 

Yet another rapidly developing related area is that of 
quantum information processing with trapped ions. Af- 



ter the first works, in which the 3- and 4-ion GH Z-state 
prj| . and 3- ion W- and GH Z-state [H| have been gen- 
erated [l2j , in recent experiments the tomography of 6-, 



13j , and the 



7-, and 8-ion W-states has been performed 
6-ion GH Z-state has been generated (l4| . 

The problem of characterization of the generated forms 
of multipartite entanglement [Hj], or more generally, of 
characterization of many-body quantum correlations is 
thus of essential importance for the investigations of 
such mesoscopic systems. One of the possible ways to 
achieve it, is to measure the total spin (or pseudo-spin) 
of atoms (or ions) and its quantum fluctuations, which 
is, of course, possible by performing state tomography. 
The central role in this approach, applied to atomic en- 
sembles, has been played so far by the, so called, spin 
squeezing parameter £ 2 , introduced by M. Kitagawa and 
M. Ueda in Ref. As it was shown in Refs. pj], El 

it provides a sufficient entanglement criterion for atomic 
ensembles. On top of that, £ 2 is particularly appreciated 
by experimentalists since: i) it has a clear physical mean- 
ing, ii) it can be relatively easy measured, iii) it is defined 
by a simple operational expression, iv) it provides a fig- 
ure of merit for atomic clocks. However, until our recent 
Letter [l9( no further investigations to relate £ 2 to other 
concepts of quantum information have been carried out. 

The present work is a substantially extended version 
of Ref. [I|| . Apart from the expanded theoretical analy- 
sis, we present here a detailed description of the ion-trap 
experiment of Ref. [13], to which output we apply our 
inequalities. 

Let us first recall that in Ref. [Il we have generalized 
and connected the concept of spin squeezing parameters 
to the theory of entanglement witnesses [2fJ, i. e. such 
observables W that have non-negative averages for all 
separable states and there exists an entangled state g 
such that tr(pVV) < 0. In order to derive the generalized 
spin squeezing inequalities, we have proposed a general 
method of expressing state averages of the appropriate 
entanglement witnesses in terms of the macroscopic spin 
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operators: 

N 1 

Jl = Y,2<> *' = 1 > 2 >3 (!) 

a=l 

(we work in the units h = 1). Here by a 1 we denote 
Pauli matrices, indices a, b, c . . . enumerate the particles 
of the ensemble, and a l a — a 1 (g> li...a...iv (hat over the 
index denotes that it is omitted). It is worth recalling at 
this place that in the standard terminology [IB] a state 
of a spin- J system is called spin squeezed if there exists a 
direction n, orthogonal to the mean spin (J), such that: 

e = 2<AJ 2 )/J < 1, (2) 

where J n = n • J. 

Our method works as follows: we begin with consider- 
ing symmetric states of N qubits first, i.e. states g that 
satisfy: 

PqP = q, (3) 

where P is an orthogonal projector onto the symmetrized 
product of individual qubit spaces Ti s = Sym(C 2 ® • • • (8 
C 2 ) (Sym denotes symmetrization) . It is known that 
for symmetric states of two and three qubits the neces- 
sary and sufficient condition for separability of a quan- 
tum state is equivalent to the positivity of partial trans- 
pose (PPT) of the state. For two qubits the PPT con- 
dition is in fact the necessary and sufficient condition 
for separability of arbitrary (also non-symmetric) states 
[2l| ; for symmetric states of three qubits this result has 
been shown in Ref . [22] . The knowledge of the necessary 
and sufficient separability criterion allows us to derive the 
complete families of generalized spin squeezing inequal- 
ities, which provide necessary and sufficient conditions 
for genuine 2-, or 3- qubit entanglement for symmetric 
states. 

Our inequalities at the same time provide a sufficient 
condition for entanglement of general, i.e. not necessarily 
symmetric, states of N qubits (23|. The results of Ref. 
[191 ] imply also that, if we somewhat broaden the stan- 
dard notion of spin squeezing @, then for spin- J systems 
represented as a collection of 2 J qubits, spin squeezing 
becomes equivalent to the bipartite entanglement among 
the qubits (see also [181 ] where the implication in one 
direction was obtained). We also derive and discuss im- 
proved w.r.t. Ref. [IH versions of somewhat simpler spin 
squeezing inequalities that provide sufficient conditions 
for genuine 3-qubit entanglement. 

To prepare the necessary data for the analysis of the 
output of the experiment from Ref. (l3| . as well as to 
show how to obtain our inequalities for concrete pur- 
poses, we present in this paper a very explicit derivation 
of the inequalities for the, so called, Dicke states [2~4| . 
sometimes also called generalized IV-states. We show 
step-by-step how to derive the inequalities probing gen- 
uine 2- and 3-qubit entanglement of this states. We also 



calculate all the necessary data data for checking 7- and 
8-qubit W states, which are of particular interest for us. 

In the part of our work dedicated to the experiment, 
we present a detailed description and analysis of the ex- 
perimental production and state tomography of 6-, 7-, 
and 8-particle IV-states of trapped ions, first reported in 
Ref. Here, we describe the details of the produc- 

tion of the states in an ion trap, dedicated to quantum 
information processing [25[ • We explain the step-by-step 
generation algorithm, that was implemented in the ex- 
periment. We then explain how the state tomography 
was performed, and show the full reconstructed density 
matrix of the 7-qubit TV-state. We analyze the experi- 
mental imperfections as well. Finally, we apply our spin 
squeezing inequalities to the experimental data to con- 
firm the presence of 2- and 3-qubit entanglement in the 
generated states. 

Let us stress that all of the proposed novel inequalities, 
analogously as the previously known squeezing parame- 
ter, i) have a clear physical meaning in terms of gener- 
alized squeezing and entanglement conditions, ii) can be 
relatively easy measured, and iii) are given by complex, 
but elementary expressions. Although in this paper we 
apply our theoretical tools to a fully restored density ma- 
trix from Ref. [131 ] . it is very important to understand 
that these tools require measurements of low order mo- 
ments of the total spin fluctuations only. Hence, checking 
of our inequalities can be relatively directly performed in 
large systems, such as atomic ensembles, where, in gen- 
eral, quantum tomography is not feasible. 

We also note that recently G. Toth has also derived 
various types of entanglement criteria based on entan- 
glement witnesses and on the uncertainty of collective 
observables, such as the total spin or energy [2(|. These 
criteria are useful to detect the, so-called, cluster states 
and many body singlet states, but they may also be used 
to detect the Dicke states discussed in this paper. 

The work is organized as follows: in Sections II and III 
we revise the derivation of 2- and 3-qubit entanglement 
criteria. We give there some calculational details as well 
as correct versions of the formulas (17), (18) and (19) 
from Ref. [191 ] . These results are very general and ap- 
ply to any system of qubits in any state: from few ions, 
through atomic ensembles to ultracold spinor gases. In 
Section IV we specify our inequalities to a concrete ex- 
ample of experimentally accessible Dicke states, and show 
how to construct our criteria in this case. In particular, 
we provide here explicit data for the case of 7- and 8- 
qubit IV-states, preparing the input data for the analysis 
of the experiment of Ref. [l3(. Section V is devoted 
to a detailed description of this experiment. We apply 
here our inequalities to the analysis of the output of this 
experiment, confirming the presence of 2- and 3-qubit en- 
tanglement. Section VI is dedicated to the construction 
of simplified witnesses detecting genuine 3-qubit entan- 
glement, impro ving the similar witnesses we constructed 
earlier in Ref. [19]. The simplified witnesses, unfortu- 
nately, are not very useful for W-states, since they detect 
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entanglement only for low number of qubits. That is why 
we do not use them for analysis of the experimental data 
of Ref . [HI . We summarize our results in Section VII. 



II. DETECTION OF BIPARTITE 
ENTANGLEMENT 

In this Section we present a detailed derivation of gen- 
eralized squeezing inequalities that detect 2-qubit entan- 
glement. The aim is to use the quantum fluctuations 
of the total spin, whose low moments can be relatively 
easily measured, as an indicator of entanglement. The 
results obtained in this Section, as well as the method 
itself, have a very general character, and can be used for 
arbitrary systems of qubits in any quantum state. 

Let us recall that a multiqubit state g possesses 2- 
qubit entanglement if for some qubits a and b the reduced 
density matrix: 



Qab = tr 



l..a..b..N 



Q 



(4) 



is entangled. The PPT criterion [2l| implies that g a b is 
entangled iff there exists a vector \tp) such that: 



tr. 



AeabH)(^\ Tl ) <o, 



(5) 



where transpose is defined with respect to (w.r.t.) the 
standard basis |0), As ip we can take any eigenvector 
of g^l corresponding to a negative eigenvalue. 

According to our general strategy we first consider 
symmetric states, as then we can obtain a convenient 
parametrization of \ip). In the 2-qubit case we can take 
advantage of the low dimensionality and use the explicit 
form of g^l . Let us first fix the basis of each qubit space 
by a z \0) = |0), a z \l) = -|1>. Then we have that: 



eo 5 S* T 
5* £l w* 

S vj ei <r 



(6) 



where e ,ei,£2,T G 
the type: 



I. It is easy to check that vectors of 



|V)=t?|00)+/3|01)+/3*|10)+7|11), 77,7 eM (7) 

are preserved by g Tl , and, since they have three inde- 
pendent parameters (we take them to be normalized, al- 
though it is not important for the condition (J5j) ) , it is 
possible to find a solution of the eigenvalue equation. 
Hence, the negative eigenvalue vector in the inequality 
([5]) must be of this form. From Eq. ([7]) it follows that 
the matrix [ip] of coefficients of \ip) is hermitian: 



(3* 7 



(8) 



and hence we can diagonalize it by some U <E SU(2) 
(modulo U(l) phase rotation): 



[ip] = ffl&U. 



(9) 



Note that due to the normalization of IV 7 ) , the eigenvalue 
matrix A can be put in the following form: 



A = 









±cos^ 



— 7T < Ifi < TT. 



(10) 



Rewriting Eq. §E§ explicitly in the basis, and using Eq. 
(|10p . we finally obtain the following parametrization from 
Ref. 0: 

|V) = U* ® C/|Vo) , IV'o) =sin||00) + cos||ll), (11) 

where U = U T , and we have fixed the overall phase. The 
parameters rj, /3, 7 from the decomposition ([7]) are now 
encoded into ip and U. Using the above parametrization 
inequality takes the following form: 

tr ab ( gab U ®U\ip )(i> \ Tl tf <0. (12) 

In order to rewrite the condition (fT2|) with the total 
spin operators ((T|), we first recall that |V'o)('0o| Tl can be 
decomposed into Pauli matrices, as it was done in Ref. 
[l9l |. Then, the adjoint action of SU(2) in the inequality 
(fT2"|) induces a 5*0(3) rotation R of the Pauli matrices: 
Ucr l U^ — K 1 (here and throughout the work we sum 
over the repeated indices). We will denote the axes of 
the rotated frame by k, 1, n. 

Since in the symmetric case we currently consider all 
the reductions g a b are of the same form, we can sum 
the inequalities (fT2|) over all pairs of qubits: J2( a b) = 

X2o=i 1 X/6=a+i> without affecting the inequality sign. 

However, before we do so, we extend \ip){'4>\ Tl from the 

space of the qubits ab to the full Hilbert space of TV qubits 

by: \ip)ab{i>\ Tl = |V)(V| Tl ® K.&..b..N- Then ' we obtain 
that: 



J2 tr AeabH)(^\ Tl ) =tr(g 



(ab) 



Ei 

(ab) 



ab{ 



(13) 



Now we can plug the Pauli matrix decomposition of 
\ip)(ip\ Tl into Eq. Q3J), and, using the identity: 



^<®^=2(J 

(ab) 



i\2 



N 

y 



(14) 



obtain the desired form of the condition ([5]), i.e. a sym- 
metric state g possesses bipartite entanglement iff there 
exist -7r < if < 7r and U G SU(2)/U(1), such that the 
following inequality holds: 



swip 



[Jl) 



(j, 2 ) 



(Jl) 



N(N - 2) 



N 
~2 

<0, 



(N - l)cos<p( J a ) 



(15) 



where all the averages are taken w.r.t. the full state g. 

In case of a general, i.e. not necessarily symmetric, 
state g observe that, if there exist — ir < ip < n and 
U G SU (2) /U(l) the same for all pairs of qubits, and such 
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that the sum (fl~3f is negative, then there must be at least 
one pair ab for which tr a ;,(p a ; ) |'0)('0| Tl ) < 0: an d hence 
the state g possesses bipartite entanglement. Thus, the 
condition (|15[) is also a sufficient condition for bipartite 
entanglement for general states. 

For a given negative eigenvalue vector \ijj) the left 
hand side of the inequality (|15[) is completely determined. 
However, we can also treat it as a function of the parame- 
ters of \if>), and as such it can be optimized. In particular, 
keeping the frame k, 1, n fixed, we can search for the min- 
imum w.r.t. ip. Let us call this minimum ip . Clearly, if 
the inequality (fT5|) is satisfied for some <p, then it will be 
also satisfied for ipo, and vice versa. Hence, it is enough 
to check the condition (fT5|) only for tpo. Performing the 
minimization, we obtain that: 



smip = 



COSi^o 



(J) 



(J?) f 



[{Jy 



(J?) 



f] 



+ (N-i) 2 (j n y< 



(16) 



(N-l)(J a ) 



\/[(JZ)+(J?)-f} 2 + ( N - 1 ) 2 ( J n) 

and the inequality (fl~5)) becomes: 



=, (17) 



N(N — 2) 



< 



N 
~2 



(Jl) 



(Ji 



N(N + 2) 



the criterion (fl8| can be simplified to: 



4(AJ 2 
N 



< 1 



4(J n 



N 2 



(19) 



(20) 



The relation of the criterion (jT5J) to the standard spin 
squeezing condition @ is the following. Spin-J state 
can be equivalently represented as a symmetric state of 
N = 2J qubits. Intuitively, spin squeezing should refer 
to the existence of non-classical correlations among the 
qubits [16]. Indeed, the criterion (|20p provides a rigorous 
proof for this intuitive picture, as, on one hand, if the 
condition {2} is satisfied, then the inequality (|20|) is sat- 
isfied as well, since in this particular case (J n ) = and 
J = N/2. Hence, spin-J squeezed states possess 2-qubit 
entanglement On the other hand, if we broaden the 
standard definition of spin squeezing and allow the 
direction n to be arbitrary, then we also obtain the con- 
verse statement: the condition (f2"0"l) implies the existence 



of a spin component </„, such that (AJ 2 ) < J/2. Note 
however, that from the condition (|2H]l it docs not follow 
that the direction of squeezing n is orthogonal to (J). 
Thus, we obtain a more general type of squeezing. In 
Section IV we will show somewhat extreme examples of 
state, for which n is actually parallel to the mean spin. 



III. DETECTION OF TRIPARTITE 
ENTANGLEMENT 

In the previous Section we have shown that the pres- 
ence of bipartite entanglement is detected by the second 
order moments of the total spin. Given our methods, it is 
natural to expect that the tripartite entanglement should 
be detectable by third order moments of J. In this Sec- 
tions we derive the corresponding generalized squeezing 
inequalities in the most generic form, valid for arbitrary 
quantum states of systems of qubits. 

As in the previous Section, we begin with consider- 
ing symmetric states first. Recall that the PPT criterion 
still works for the tripartite reductions g a b c of such states 
[12], and there are two families of potential negative- 
eigenvalue- vectors of g T a \ c [TJ|, H3 : 



+ (N - l) 2 (J n > 2 . 
(18) 

As a result, we arrive at the |19j : 

Criterion for bipartite entanglement. If there ex- 
ist mutually orthogonal directions k, 1, n such that the 
inequality H8\) holds, then the state g possesses bipartite 
entanglement. For symmetric states the above condition 
is both necessary and sufficient. 

In the latter case, due to the equality: 



|V.'> = A®B(g>B\GHZ 3 } 
|V> = A®U ®U\W 3 ) . 



(21) 
(22) 



Here, matrices A, B G SL{2,C), U G SU{2), and 
\GHZ 3 ) = (1/V2)(|000) + |111», \W 3 ) = (1/V3)(|001) + 
|010) + I LOO)). The action of SX(2,C) on the Pauli ma- 
trices in the decomposition of (i[>\ Tl now induces re- 
stricted, i.e. orientation and time-orientation preserving, 
Lorentz transformations: 

A*a^A T = A" v a v , Ba^B* = L^ v a u , a = 1 (23) 

(Greek indices run through ... 3), and the PPT condi- 
tion takes the following form: 

tr afcc ( g a bc | ip)(tp | Tl ) = \K aPl ii abc (g abc a a <g> a p <g> a 1 ) < 
(note the summation convention). Tensor K a pj reads: 

K aPl (k 1 L) = k\L%L\+k\L%L\+k l a L 1 p L\ 
A 1 a L 2 L\ + 2A 2 a L\ L\ ) (25) 



-2A d L 



for the G-ffZ-family ((21]) . and [2j: 

Kapyi-A, R) = ||3A° Q i? /3 i?° 7 — 3A 3 a R 3 pR 3 J 
-2A° n R° m R 3 ^+A 3 a R° pR° y — A° ct i? 3 /3 i? 3 7 (26) 



~2A 3 Q i? (/3 _R 3 7) +4A 1 Q i? (/3 i? 1 7) +4A 1 Q _R 1 ((3 i? 3 7) 
-4A a R ^R ^ - 4A 2 Q i? 2 (/3 i? 3 7) + 2A° a i? 1 /3 i? 1 7 

\-2A 3 a R 1 ^R 1 y + 2A 3 a R 2 pR 2 7 + 2A° Q i? 2 /3 i? 2 7 | 



for the W-family (|22)l . Here, i? M „ is the four-dimensional 
embedding of the rotation generated by U from Eq. 
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and the round brackets ( ) around indices denote sym- 
metrization, e.g. = (A ij + A^ l )/2. Note that the 

relativistic notation is used only for our convenience. We 
could have as well put all the indices at the same level as 
we are not going to lower or rise them with the Minkowski 
metric. 

Next, we sum the inequalities over all triples of 
qubits: £ (afcc) = Y,a=i : E^+iEfLfe+n just like we 
summed the inequalities (|12|) in the previous Section: 



^ K aj3l tT: abc {Q abc a a (g) a 13 <g> a 1 ) = 

(abc) 

^W^S ?® *®^)- (27) 

{abc) 

Because of the symmetry condition @, we can rewrite 
Eq. (J27J) as follows: 



K a(3l tY(gP ^ o% ® erf ® ^P 1 ), 



(28) 



(abc) 



and observe that due to the action of P, we can substitute 

a a® a b® a c w i tn the symmetrized product a a a ®(7 b ®crl) ■ 
This finally allows us to rewrite Eq. (|2"8"]) with the total 
spin operators J 4 (supplemented by an artificial "time- 
component" J° = (N/2)l for compactness of the nota- 
tion), because of the identity: 



(29) 



The symmetrization above is taken w.r.t. a/?7 and 7 /it 
separately and square brackets [ ] around Greek indices 
denote antisymmetrization, e.g. A^ = (A tj — A jl )/2. 
The constants f a ^ are defined through: cr^o" — /^er 7 . 



Their numerical values are as follows: f 0< % — f a0 a 



3i 



f \ = iEi^'^ + ^^a- Substituting Eq. ((29j> into 
Eq. (EHD leads us to the 

Criterion for tripartite entanglement. A symmet- 
ric state g possesses genuine tripartite entanglement iff 
there exist two restricted Lorentz transformations A, L, 
or a restricted Lorentz transformation A and a rotation 
R, such that: 



X{g) ee K {aPl) [2(rjP.r)~Zt^(J^J^) (30) 



< 



holds, with K a pj given by Eq. i25\) . or by Eq. 

For a general state g we could, as in the previous Sec- 
tion, generate a sufficient entanglement condition by ap- 
plying the same witness \ip}(ip\ Tl , with \ip) given by Eq. 
(|2"T|) or by Eq. (f2"2"j) . to all tripartite reductions g abc . 
However, then we cannot use the symmetry arguments 
like we used in Eq. (|2"8"]) . and directly apply the identity 



(129)) . Instead, we construct from the families (|2ip and 
(|2"2"| different witnesses, given for the W^-family (f2"2")l by: 

i{ (AUU\W 3 )(W 3 \AWtf) Tl 
+ (UAU\W 3 )(W 3 \U^A^U^ T2 
+ (UUA\W 3 )(W 3 \U^UU^ T3 Y (31) 

(we have omitted tensor product signs, (g>, for compact- 
ness), and analogously for the GHZ- family (f2"Tj) . We 
then apply the witnesses (f3Tj) to all tripartite reductions 
of g, which effectively leads to the substitution of K a ^ 
by K [aPl) in Eq. $27} [HJ. Hence, we can use Eq. 1(25]) 
again and arrive at the condition (|30p . 

The price to pay, apart from the mere sufficiency of 
the condition (|30p , is that the witnesses (f3"Tj) make no dis- 
tinction between biseparable tripartite reductions (which 
now are not forbidden by the symmetry) and genuine 3- 
qubit entangled ones, and hence, the inequality (|30)) indi- 
cates only general 3-qubit entanglement. However, note 
that the set of all biseparable states is closed, and hence 
each genuine 3-qubit entangled state possesses an open 
neighborhood consisting of only genuine 3-qubit entan- 
gled states. Thus, the criterion (|3T))) also detects genuine 
3-qubit entangled states in some open vicinity of sym- 
metric states, but the size of this vicinity is a priori not 
known (the same remark applies to the criterion (|2H)) as 
well) . We will partially solve this drawback using another 
witnesses in Section VI. 



IV. AN EXAMPLE - DICKE STATES 

There exists a famous and experimentally accessible 
family of symmetric states, called Dicke states [24j : 



I* 



N,k) 



1 000 ... 0) + perm) , (32) 



('perm' stands for all possible remaining permutations), 
which are generalizations of iV-qubit VF-states |H^v): 



W N ) = |tf 



NA) 



(33) 



Dicke states correspond (in spin 1/2 language) to a fully 
symmetric flip of k out of N spins. Such states appear in 
many physical processes, such as superradiance, super- 
fluorescence. As we already mentioned, they can also be 
realized with photons (where the qubits are encoded in 
the polarization degree of freedom), or trapped ion sys- 
tems (where the qubits correspond to two internal states 
of ions). In this Section we specify and apply our general 
results of the previous Sections to the Dicke states. We 
explicitly construct for \^N,k) the inequalities (fT8|) and 
(I30p . In particular, we derive all necessary expressions for 
the analysis of the experimental data on 7- and 8-qubit 
VF-states, which we will perform in the next Section. 



For practical reasons, we choose the number of excited 
qubits k to be smaller than the integer part of N/2. Also 
note that alternatively the states (|32|) can be defined as 
the eigenstates of the total angular momentum: 



states: 



\* Nfk ) = \N/2,N/2-k). 



(34) 



We first consider 2-qubit entanglement. The reduced 
2-qubit density matrices all have the following form: 



Q2 = 

where 



-i 

(c o |OO)(OO| + Cl |ll)(ll| + 2 C+ |0 + )(^ 



|^+) = ^(|01) + |10)), 



(35) 



(36) 



and the coefficients are given by the following binomials: 



Co 



N-2 

k 



Cl 



N-2 
k-2 



N-2 
k-l 



In the basis |00), |11), |01), 1 10} the partially transposed 
matrix g^ 1 is given by: 



/ co c+ 

c+ ci 

c+ 

V c+ 



(38) 



In the generic case, when all the constants from Eqs. (13"?| 
are non-zero, g^ 1 has one negative eigenvalue: 



as coCi — c? + < 0, and hence the states ([3"2l possess bi- 
partite entanglement. The normalized eigenvector corre- 
sponding to A_ is given by: 



1 



cp - Cl 
2c+ 



(|00)-i|ll», 



cp - Cl 
2c+ 




(40) 
(41) 



We see that \ip) is already in the Schmidt decomposed 
form w.r.t. the chosen basis, and hence no unitary rota- 
tion U is needed. As that rotation was the only ingre- 
dient needed to construct the spin squeezing inequalities 
(fl"8)) and (|20|) (because the angle ip is minimized over), 
we simply put k, 1, n = x, y, z in them. 

Although in theory both inequalities (fT8"|) and P0|) are 
equivalent, and we could use the latter due to simplicity, 
the inequality to be measured is rather (fT5|) as in real- 
life experiments one does not obtain perfectly symmetric 
states. Using Eq. (f3"4")) we find that for the perfect Dicke 



N{N -2) _ N(N-l) 



-Nk + k 2 , (42) 



((^) 2 > + ((^) 2 >"f n2 



+ (N — 1) 2 {J Z Y 



{Nk-k 2 ) 2 + 



(N -l) 2 (N -2k) 2 



(43) 



For the experimentally interesting examples of the 7- and 
8-qubit W-states \Wr), \Wa), the expressions (1421) and 
([4*3")) take the following values: 15.000 and 16.155 respec- 
tively for |W 7 ); 21.000 and 22.136 respectively for |W 8 ). 

Let us now proceed with the analysis of tripartite en- 
tanglement. All tripartite reductions are of the form: 



£':s 



N\~7 

J (k |000)<000| +Ki|lll)(lll| + 
+3uj\W 3 )(W 3 \ + 3u;'\Wi)(Wi\) , (44) 



(45) 
(46) 



where \W^) = 1/V3 (|011> + |101> + |110)), and 



K 



N -3 
k 

N -3 
k-l 



N -3 
k-3 
N -3 
k-2 



In the basis |000), |110), |101), |010),|001), |111),|100), |011) 
the partially transposed matrix g^ 1 reads: 



A- = \ (J) (co + ci - .^(co-c^+Acl) , (39) ^ _ /jyx 



LU to' iO 1 








V o 



\ 





to lu J 

ui u lu' 

ui' J ki 

w 

lu' J 



(47) 



In the generic case it has two negative eigenvalues: 

M- = ^(T) («o + 2w ' - V(«o - 2^') 2 + 8^ 2 ) , (48) 
- = \(je) ( Kl + 2w — -\/( K i ~ 2uj) 2 + 8lu' 2 ) (49) 



(because kquj' < ui 2 and K\U) < to' 2 ), and thus the states 
\^N,k) possess tripartite entanglement as well. Since 
there are two generically different negative eigenvalues, 
there will be two different spin squeezing inequalities 
([30)) . As before, we will generate them from the eigen- 
vectors corresponding to fi- and //_, which read: 

\ip) = |000)-a|l)®(|01) + |10)), (50) 
|V/} = |lll)-a'|0)® (|01) + |10)), (51) 
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where: 



o 



k — 2lu' 

Ki — 2lu 
Auji ' 



I/k - 2uj'\- 
V V — 4^ — J 



/Ki - 2ui\ 2 
V 4w' ) 



1 

2' 



(52) 
(53) 



The vectors ([50)1 , f|5 1 [) are not normalized, as the norm 
is irrelevant for the PPT condition (J2H). After proper 
rescaling, \tp} and \ip') can be rewritten in the desired 
form d22|): 



|V;) = A® 1® 1|W 3 ), 



where A, A' S SL(2,C) are defined as follows: 
A- ± ( 
A' = ± 



4= 

-y/a ' 
iVa 7 







(54) 
(55) 



(56) 
(57) 



Before we proceed with the construction of the inequal- 
ities pop , let us note that having the explicit forms of the 
negative eigenvalues and the corresponding eigenvectors 
of Q3 1 it is straightforward to calculate the sum over all 
triples of qubits (j2"4")l . It is just given by: 



E(abc) ^abc(QabcW)(lp\ 



IV 



N 



N\ 2c? 



3o 



(58) 



for and |?/>), and by the analogous expression for fjf_ 
and |V). However, our goal here is to express Eq. 
using total angular momentum, in order to make it exper- 
imentally available and connect it with the spin squeez- 
ing. 

Hence, following the procedure described in Section 
III, we first have to find the Lorentz transformations and 
rotations generated by matrices from Eqs. (|54|) and (|55|) . 
These transformations are the following: matrix (|56p gen- 
erates, according to Eq. (j2"3")l , the rotation by tt around 
?/-axis, followed by a boost along z-axis: 



A(A) 

P 



( 7 -7/3 
0-100 
10 
V 7/3 -7 



a 2 -l 



1 



1 



, 7 = 



^1 



(59) 
(60) 



Obviously the identity operator 1 from Eq. (|54[) gener- 
ates the trivial rotation, so we have in this case R = 1. 
Matrix (|57[) generates the rotation by n around z-axis, 





FIG. 1: (Color online) The (interpolated) plots of the l.h.s. 
of the inequality (|30p corresponding to the eigenvectors \ip) 
(left) and (right). 



followed by a boost along it: 
A(A') 







( i iff 
0-100 
00-10 
7' 

' 2 - 1 , 1 



' 2 + l' 7 y/T^W 



(61) 
(62) 



while (t x from Eq. ([55)1 generates the rotations by it 
around a;- axis: 



R(a x ) 



/ 1 
10 
0-10 

V -1 



(63) 



in the spaces of the second and the third qubit. 

Next, from the matrices A(A), R = 1, and A(A'), 
R(a x ), we construct two copies of the tensor K a p^ : ac- 
cording to Eq. (|26p . Finally, having K a ^, we construct 
the corresponding parameters X(Wjv.fc), defined in Eq. 
(|3TJ|) . and check the 3-qubit spin squeezing inequalities. 
The resulting expressions are lengthy but straightfor- 
ward, and hence we will omit them here. Let us stress 
that for the ideal, generic Dicke states we obtain two 
independent inequalities, and both of them must be sat- 
isfied. Fig. [Qshows the plots of X(^> M,k) as a function of 
N and k [3lf. 

Let us now analyze the ./V-qubit VF-states \Wn) of Eq. 
(f3"3"| . In this case, from Eqs. (|4"5| and (|4"6"|) we see that 
Kq = N — 3, w = 1, and K\ = u>' = 0. Substituting 
this constants into Eqs. (|48[) and ([49]) . we obtain that 
there remains only one negative eigenvalue of g^ 1 given 
by /!_. As a consequence, states \Wm) lead to only one 
spin squeezing inequality, generated by the matrix A{A) 
from Eq. (f59|) and the trivial rotation R = 1. The pa- 
rameter a from Eq. (160|) is now equal to: 



iV-3 



iV-3\2 



(64) 



8 





(c) 

|0>J0> l '>J 0> 



l 0> m 



M>J 1 > 



FIG. 3: (a) Level scheme of Ca + . (b) Schematics of the 
two lowest levels of the harmonic oscillator describing the bus 
mode, (c) Joint energy level diagram of the electronic qubit 
levels {|1),|0)} and the phonon numbers of the ion's motional 
mode used for entanglement generation {|0} m ,|l) m }. Carrier 
transitions are marked as solid arrows, the blue sideband tran- 
sition as a dashed arrow. Note that the |0} m |0}-level does not 
couple to the blue sideband. 



FIG. 2: The plot of the parameter X as a function of the 
amount of noise for noisy W^-states of N = 7 qubits (solid 
line) and N — 8 qubits (dashed line). 



For the state \Wr) we obtain from the corresponding 
formulas that: 



A(A) 



1.337 -O.l 
0-10 
10 

0.888 -1.337 



(65) 



and X{W 7 ) = -44.04. 

For the state |Ws), the corresponding matrix is given 
by: 



A(A) 



( 1.529 



V 1.157 




(66) 



and the parameter X(Ws) = —59.88. 

To better understand the meaning of the above values 
of the parameter X^n^), let us briefly consider a less 
idealized situation and mix the states (1321) with the white 



noise: 



Q = p\^N,k)(^N,k\ + ( 1 -P)^W 



(67) 



We then calculate the parameter X(g) as if the state (j6"?| 
were an experimental output: we calculate the averages 
of the spin operators in Eq. (|30[) using the density ma- 
trix (|67| . while plugging the tensor K a ^ calculated for 
the ideal Dicke states. Thus, X(g) = pX(^ N , k ) + (1 - 
p)X(l/2 N ). The results for the states \W 7 ) and |W 8 ) 
are presented in Fig. [21 



V. EXPERIMENTAL RESULTS 

The aim of this Section is to apply the tools developed 
in the previous Section to the recent experiment of Ref. 



[13] . In this experiment, 7- and 8-qubit M^-states have 
been produced in an ion trap, dedicated to quantum in- 
formation processing [25| . We begin this Sections by pre- 
senting necessary details of the experiment, and follow by 
applications of our generalized squeezing inequalities. 



A. Description of the experiment 

Strings of up to eight 40 Ca + ions are held in a lin- 
ear ion trap capable of storing the ions for several days, 
a time sufficiently long for creating an entangled state 
more than 10 6 times. The qubits are encoded in super- 
positions of the S1/2 ground state and the metastable 
D 5 / 2 state of the Ca + ions (lifetime of the -D5/2 level: 
t 1.16 s). For the atomic level scheme, we refer to 
Fig. [3^. Each ion in the linear string is individually ad- 
dressed by a series of tightly focused laser pulses on the 



|1> = 5 1/2 (mj 



-1/2) 



|0> = As/2 ( m 



-1/2) 



quadrupole transition with narrowband laser radiation 
near 729 nm. Depending on its frequency, the laser cou- 
ples either the states |rt) m |l) <-> |ri) m |0) (carrier pulse) 
or the states |n) m |l) |n+l) m |0) (blue sideband pulse, 
laser detuned by +uj z w.r.t. the atomic transition, see 
Fig- (St)- Here, n denotes the vibrational quantum num- 
ber of the ion string's center-of-mass motion. Via side- 
band cooling and optical pumping, the ions are prepared 
in the |0) m |ll • • • l)-state. 
The TV-ion VF-states: 



\0} m \W N ) 



1 



£lx<>> 



\Xi) 



N ■ 

\0} m \x N ...Xi) 

1, if k = i 
0, if k ^ i 



(68) 



(note the reverse ordering of the qubits) are created by 
applying the sequence of laser pulses shown in Table [I] 
to the ions. First, the |0) m |lll • • ■ 1)— state is prepared 
by N 7r-pulses on the carrier transition applied to ions 
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TABLE I: Creation of a |Wjv)-state (N = {6,7,8}). The numbers within the state vector refer to the phonon excitations of 
the center-of-mass mode of the ion crystal. The electronic states are labeled by |1) and |0). Rn{6) denotes a carrier pulse 
of length 6 applied to the ion n, Rt(9) a blue sideband pulse, (il) • • • (i3) mark initialization steps, (1) ■ ■ ■ (N) the actual 
entangling steps. Note that we count the atoms from right to left. 



(il) 
(i2) 

(i3) 
(1) 
(2) 



(N) 



|0) m |lll---l) 



|0) m |000---0) 

Check state via fluorescence 



|0) m |000--0> 

Check state via fluorescence 



^|0> m 100 ■■■0) 

i?t(2 arccos(l/\/jV) 



^=|0> m 100 ■■■0> + 

H+ _ 1 (2 arcsin(l/VJV-l) 



V^V 



U0) m 100 •■•()> + -£ 



|l) m |000---0) 



|010 • ■ • 0) + 



1 

H+(2 arcsintl/v^T) 



100 • ■ ■ 0) + 



010---0) + • 

^|0} m 100---0> + ^|0> m |010---0) + . 



|l) m |000---0) 



Vff 1 



+ _^|l) m |000---0) 
+ -*=|0> m |000---l) 



#1 to #N. Then, laser light coupling the |1) state reso- 
nantly to the short-lived excited state P\/2 projects the 
ion string on the measurement basis. Absence of fluo- 
rescence reveals whether all ions were prepared in |0). 
Similarly, we test the motional state with a single blue 7r 
pulse. Absence of fluorescence during a subsequent de- 
tection period indicates ground state occupation. This 
initialization procedure can be viewed as a generalized 
optical pumping with the target state |0) m |ll ■•■!). If 
both checks were successful (total success rate > 0.7), 
we continue with the | W)-preparation at step (i3) in 
Tab. [J to create the state |0) m |10- • -0). The entangling 
procedure starts by moving most of the population to 
the 1 1) m. 1 000 - • • 0) with a blue sideband pulse of pulse 
area On — 2 arccos(l/\/^V) leaving 1/N of the popula- 
tion back in |0) m |100 • • • 0). Now, VF-states arc efficiently 
generated by redistributing the |l) m |0 ... 0) state popu- 
lation equally among the states |0) m |0 . . . 01^0 . . . 0),i = 
1 . . . N— 1. This is achieved by N — 1 blue sideband pulses 
of pulse length 0i = 2 arcsin(l/\/Ar — i). Note that for 
an ion string in the motional ground state, blue-sideband 
pulses acting on an ion in the |0)-state have no effect. We 
note that this production method scales quite advanta- 
geous, as the required sideband pulse area increases only 
logarithmically. Therefore, even large W^-states can be 
created quite efficiently. 

Furthermore, the space spanned by the states \xi) from 
Eq. (|68[) is decoherence-free w.r.t. the collective de- 



phasing, which is the main decoherence mechanism in 
the setup [32|. Therefore, the VF-states are quite robust 
and long-lived. In addition, even during the creation of 
a W-state, the energy difference between all superposi- 
tions never exceeds more than the one of a single qubit. 
Thus, the requirements to laser frequency and magnetic 
field noise for a high-fidelity generation of W-states are 
rather modest. We discuss the imperfections in a sepa- 
rate section below. 



B. State tomography 

Information about the A^-ion quantum state is ob- 
tained by exciting the ion string on the S\/2 <-* P\/2 tran- 
sition (see Fig. [3^) , and detecting the ion's fluorescence 
spatially resolved with a CCD camera state 25 1. The 
measurement of an ion's fluorescence amounts to mea- 
suring the Pauli matrix a z , if |0) and |1) are identified 
with the eigenstates of a z . The measurement of a x (<j y ) 
is accomplished by applying a suitable 7r/2 carrier pulse 
to the ion prior to the state detection [33] ]. 

To verify the entanglement of the produced state, a 
measurement of a witness operator, yielding a negative 
expectation value, would be sufficient in principle. How- 
ever, the optimal witness is a priori not known. There- 
fore, it can be advantageous to get as much information 
as possible about the produced quantum state. Full in- 
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formation on the TV-ion entangled state is obtained via 
quantum state reconstruction. For this we expand the 
density matrix in a basis of observables, and measure 
the corresponding expectation values. For the basis, we 
choose tensor products of Pauli matrices: a]^ ® . . . ® o^ 1 
(note the reverse ordering). We use 3 N different bases 
and repeat the experiment 100 times for each basis. For 
N = 8, we need thus 656 100 experiments and a to- 
tal measurement time of 10 hours. We follow the iter- 
ative procedure outlined in Ref. [34| for performing a 
maximum-likelihood estimation of p. Other reconstruc- 
tion methods would also be possible 35| . The procedure 
ensures also positivity of the reconstructed matrix. The 
resulting matrix for the state \Wr) is displayed in Fig. [H 
the numerical values are published in the on-line material 
of Ref. pjjj. 

A Monte Carlo simulation is used to estimate uncer- 
tainties in the density matrix elements, and in quanti- 
ties derived from it, that are due to quantum noise in 
the state reconstruction measurements: starting from the 
reconstructed density matrix, we simulate the measure- 
ment process and reconstruct up to 100 times the density 
matrix from these simulated measurements. From the 
set of reconstructed density matrices, the spread in the 
expectation values of the observable of interest can be 
estimated. For density matrices close to pure states, we 
observe that the purity of the reconstructed matrices of- 
ten slightly decreases (for the \W) states by about 2%). 
Therefore, we conclude that the reconstruction process 
rather underestimates the entanglement in the experi- 
mentally produced quantum states. 



C. Experimental imperfections 

For an investigation of the experimental imperfections 
and scalability, we simulate the preparation procedure 
by solving the Schrodinger equation with the relevant 
imperfections. 

The four major sources of deviations from the ideal 
W-states are addressing errors, imperfect optical pump- 
ing, non-resonant excitations, and phase noise (laser fre- 
quency and magnetic field noise). For the large W-states, 
we approximate the ions as two-level systems and include 
only the first three levels of the center-of-mass excitation. 
For a serious analysis of the imperfections this is by no 
means sufficient as e.g. no environment is included. Still 
the simulation time for the generation of a |W / s)-state 
under these idealized conditions is already 20 minutes 
on a 3 GHz processor using Matlab. As the computa- 
tional time for the simulations scales with A N , it is quite 
demanding to include a reasonable environment or even 
use a density matrix approach. 

The fidelity reduction of \Wq) for the different imper- 
fections are as follows: 0.1 (addressing error), 0.07 (off- 
resonant excitations), 0.04 [laser frequency noise (200 Hz 
rms)]. We note here that as opposed to e.g. experiments 
on teleportation [3(| or the Cirac-Zoller controlled-NOT 



|25j |. phase noise (caused by the laser frequency noise or 
magnetic field noise) contributes here much less. Another 
possible error source is imperfect ground state cooling. 
Intensity noise of the 729-laser (AJ max /7 « 0.03) does 
not contribute significantly. Finally, we experimentally 
observed non-ideal optical pumping which can result in 
a reduction of 0.02 of the fidelity per ion. For N > 6, 
we therefore minimize the errors due to optical pump- 
ing and a part of the addressing errors by checking the 
initialization procedure with a detection sequence (see 
Tab. |T|. Due to this improvements, the addressing error 
reduces the fidelity of the \ Wq) state by only 0.05 and the 
optical pumping errors are basically excluded. Further- 
more, we switched the blue-sideband pulses adiabatically 
w.r.t. the trap frequency, such that Fourier components 
at the carrier transition do not lead to off-resonant excita- 
tions. Taking this new situation into account, the fidelity 
should be of the order of 0.91. Even though it is hard to 
estimate the expected fidelity for N — 8, it seems that 
the discrepancy between the model and the experiment 
is even larger for N = 8. A small part of these discrepan- 
cies could be due to the quantum projection noise in the 
measurement process as described in the section on state 
tomography. However, looking at the density matrices in 
detail, we observe that the 1 000 • • • 0) state seems always 
quite strongly populated, especially for large N. So far 
we have no good explanation for this. 



D. Evaluation of the data 

In Ref. [l3| it has already been shown that the states 
are genuine multipartite entangled, multipartite distill- 
able and also that all the reduced two-qubit states are 
entangled. Now we want to apply our criteria to the 
experimental density matrices g ex . 

We begin with the 7-qubit states. In this case, the fi- 
delity of the produced states was F7 — 0.763. To check 
the presence of bipartite entanglement, we use the in- 
equality (fl~8|) rather than (f20|) . as the experimental states 
are not symmetric due to the experimental imperfections 
described in the previous subsection. According to the 
theoretical analysis of Section IV (c.f. formula (|4"0"]) ). the 
frame directions k, 1, n = x, y, z. We find that: 

N(N - 2) 



((J Z Y 



14.666 ±0.016. 



<(J»)3) 



N 



+ (N — l) 2 (J z ) 2 



15.148 ± 0.023. 



which clearly proves the presence of bipartite entangle- 
ment in the produced states. 

Let us move to the tripartite entanglement. We evalu- 
ate X(g ex ) using the Lorentz matrix (|65|) . We find that: 

X{gex) = -24.937 ± 0.202, (69) 

and hence the spin squeezing inequality (|30p is fulfilled. 
However, as we mentioned at the end of Section III, the 
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FIG. 4: Absolute values of the reconstructed density matrix of a |WV)-state as obtained from quantum state tomography. 
Ideally, the dark entries should all have the same height of i, the bright bars should vanish. 



validity of the inequality (|50|) only proves the presence of 
some form of tripartite entanglement and a priori we do 
not know if it is genuine 3-qubit entanglement. 

Let us now discuss the eight qubit case. Here, the 
experimentally reached fidelity was F$ = 0.7215. The 
evaluation of the bipartite criteria yields: 



((J Z Y) 



N(N - 2) 



= 20.462 ±0.007, 



<(J-) 2 ) + ((^) 2 )-y] 2 + (iV-l) 2 <J-) 2 
= 20.838 ±0.009, 
and the tripartite criterion gives: 

X{g ex ) = -29.017 ± 0.2623. (70) 
Thus both criteria detect entanglement again. 



VI. SIMPLIFIED CRITERIA FOR THE 3-QUBIT 
ENTANGLEMENT 



The general form of squeezing inequalities of Sections 
II and III is complicated and remains such even when 
specified to Dicke states. It is therefore desirable to de- 
rive alternative inequalities, which are weaker, but have 



a simple form. In fact, in Ref. [19( it was proposed to use 
less general witnesses, developed in Ref. [371] . than those 
provided by the PPT criterion for the 3-qubit case: 



W GHZ = h-\GHZ 3 )(GHZ 3 \ 
Ww x = ~l-\W a )(W 3 \, 



nv 2 



-1-\GHZ 3 )(GHZ 3 \, 



(71) 
(72) 
(73) 



where now we allow the vectors \GHZ 3 ) and \W 3 ) to be 
defined in an arbitrary frame k, 1, n, the same for all three 
qubits. Apart from the simplicity, the advantage of such 
an approach over the general criterion (j3"0|) is that the 
above witnesses detect genuine 3 qubit entanglement in 
generic states. 

We derived the spin squeezing inequalities correspond- 
ing to Wghz , Wh/j , Ww 2 using the same technique as in 
Section III: we expressed the sums 



^ tr abc (g abc Wabc) — tr vt , ^ 

(abc) (abc) 



Q Y Wabc) 



(74) 



with the total spin operators |T]). However, instead of us- 
ing the general formula (|2"9"1) , we calculated explicitly the 
occurring products of Pauli matrices (or in other words 
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we used special cases of Eq. (1291) '). This led us to the fol- 
lowing sufficient criteria for the GHZ-type entanglement 



N - 2 



i 

3^ K/ ' 2 
AT(AT- 2)(5JV- 2) 



24 



<0, 



(75) 



and for the GHZ- or W-type entanglement: 



(Jl) - 2<JiJ n Ji) - 2( J k J n J k ) 



TV- 2 



(■£» 



iV 2 - 4iV + - 



<o, 



7V(7V-2)(137V-4) 
24 

\{Jl) + (JiJkJi) - ^<^> + \w 



■(Jn) 

(76) 
(77) 



The witnesses (|71ti73p still have a disadvantage that in 
the sums Yliabc) W a bc, the identity gives the dominant 
contribution, and hence the bigger the system the less 
sensitive the witnesses become. One possible method to 
partially overcome this problem is to project the wit- 
nesses (I7H73P onto the symmetric subspace of the space 
of three qubits: 



W GHZ = jP 3 - \GHZ 3 )(GHZ 3 \ 7 (78) 
W Wl = ~P 3 -\W 3 ){W 3 \, (79) 
Ww 2 = ~P 3 - \GHZ 3 ) (GHZ 3 \ , (80) 



where: 



P 3 = |000)(000| + |111)(111| + \W 3 )(W 3 \ + 

(81) 

The factor 4/9 in the definition (f79|) is the maximum 
overlap between \W 3 ) and symmetric separable states 
(there are no symmetric biseparable states due to the 
symmetry) [38j . The criteria that such improved wit- 
nesses lead to, read respectively: 



\{Ji) + (JiJkJi) + ^(Ji + J") + \w 



<o, 



N(N -2)(N -4) 
+ 12 

- 2<J 1 J n Ji) - 2(J k J n J k ) 

7JV(JV-2)(iV-4) 



A^ 2 - 47V + , 



72 



1 



< 0, 
7V-2 



(82) 
(83) 



1, s 7V(7V- 2)(7V-4) 

iw+^r — - <0 - 



(84) 



The potential advantage of using Wghz , Wwi , WV 2 
instead of Wqhz , V^Wx , %! manifests itself only for 
non-symmetric states. For symmetric states, both fami- 
lies give the same results (apart from Wwi due to the fac- 
tor 4/9), as we can always substitute g a bc with P 3 g a bcP?, 
in Eq. (|74|) . 

Let us apply the above witnesses to the Dicke states 
I^JV.fc) of Section IV. As one can easily see from Eq. (|44|) . 
only WVi and Wwi have a chance to detect genuine tri- 
partite entanglement, but not for all N and k. For exam- 
ple, for |Wjv), WVi detects entanglement only for N < 4, 
and Wwx ~ only for A^ < 6. 



VII. CONCLUSIONS 

We develop in more detail the novel method of detect- 
ing entanglement in multiqubit systems, first introduced 
in Ref. [19( , and apply it to the output of the recent ion 
trap experiment from Ref. (l3l |. We also present some 
novel details regarding that experiment. The detection 
method is based on the use of entanglement witnesses 
together with the concept of spin squeezing. We show 
in detail and on the example, how to obtain sufficient 
(or necessary in sufficient for symmetric states) entangle- 
ment criteria for detecting 2- and 3-qubit entanglement 
(or genuine 3-qubit entanglement for symmetric states 
and states sufficiently close to them). We use them to 
analyze concrete experimental data. We also provide a 
novel, mathematically exact justification for the intuitive 
picture linking the presence of spin squeezing with the 
non-classical 2-qubit correlations in the system. There- 
fore, our criteria generalize the standard notion of spin 
squeezing as the measure of entanglement in multiqubit 
systems. 

As the concrete example we study the family of Dicke 
states and show step-by-step how our method works in 
practice. We obtain ready-to- use expressions (f6"5*)) . (|6"6"|) 
and then apply them to the study of the experimen tally 
generated 7- and 8-qubit atomic VF-states from Ref. [13| . 
In the experimental part, we explain in detail the full 
state-creation algorithm (Table [TJ) and the state tomog- 
raphy procedure used in the experiment from Ref. [13| . 
We also present the reconstructed 7-qubit W-state in Fig. 
3] and provide the analysis of the experimental imperfec- 
tions. 

Finally, we provide improved sufficient criteria — Eqs. 
(|82j|84[) . detecting genuine 3-qubit entanglement, suitable 
for macroscopic systems. 
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